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Abstract 

We study the expectation value of a Polyakov-Maldacena loop that wraps the 
thermal circle k times in strongly coupled = 4 super Yang-Mills theory. This 
is achieved by considering probe D3 and D5 brane embeddings in the dual black 
hole geometry. In contrast to multiply wound spatial Wilson loops, nontrivial 
dependence on k is captured through D5 branes. We find N^"^/^ corrections, 
reminiscent of the scaling behaviour near a Gross- Witten transition. 



1 Motivation and Introduction 



1.1 Phase structure of AT = 4 SYM theory 

The phase structure of A/" = 4 super Yang-Mills (SYM) theory at large on S*^ x 5*^ 
is of interest from both gravitational and field theoretic perspectives. The AdS/CFT 
correspondence P places both the weak and strong coupling regimes of the theory 
within computational reach, suggesting the possibility of fully mapping out the phase 
structure as a function of temperature and coupling. 

At strong coupling it was observed j21 E] that the Hawking- Page transition in ther- 
mal AdS jlj, which occurs at a critical temperature, should be interpreted as a first 
order deconfinement transition in field theory. 

At weak coupling there are also known to be phase transitions as a function of 
temperature |SIE1IZ|- In the free theory there is again a first order transition at a critical 
temperature whilst at finite small coupling there are two possibilities distinguished by 
the sign of a currently unknown coefficient in the effective Lagrangian. Either there 
is a single first order transition or, as one increases the temperature, a second order 
transition followed by a third order transition. 

The remaining question is the interpolation between weak and strong coupling 
physics. A minimal interpolation suggests that the low and high temperature regimes 
are separated by at least one phase transition at all couplings Nonrenormalisation 
theorems in the low temperature phase [Sj may suggest that the physics there depends 
smoothly on the coupling. The high temperature behaviour is less clear. In a recent 
work P we emphasised the tension between perturbative plasma physics correlators 
which generically contain many threshold branch cuts and gravitational computations 
of the same correlators at strong coupling which generically contain poles in frequency 
space. 

The order parameter used to discuss these transitions is the eigenvalue distribution 
of the Polyakov loop, the SU (N) holonomy around the Euclidean time circle 



The eigenvalue distribution may be obtained from the expectation values of the trace 
of powers of the Polyakov loop (Trf/'^). In the low temperature phase the distribution 
is uniform while at high temperatures it is not. The possible third order transition we 
mentioned is similar to the Gross- Witten transition [HI El , in which the nonuniform 
eigenvalue distribution develops a gap. 

The eigenvalue distribution of the Polyakov loop at strong coupling is not known and 
it is the purpose of this work to begin to approach the problem. An important question 



If = p ^if AodT 
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one would like to address is whether the non-uniform high temperature eigenvalue 
distribution at strong coupling is gapped or not. An ungapped distribution would 
signal a third order phase transition as a function of coupling in the high temperature 
phase 

1.2 Multiply wound Polyakov-Maldacena loops 

The AdS/CFT correspondence does not grant us access to the Polyakov loop ((H) at 
strong coupling but instead allows us to compute the Polyakov-Maldacena loop fHl 

JJ ^ p ^i§[AQ~i^ie'(T)]dT ^ j^2) 

where $/ are the six adjoint scalar fields of the J\f = 4 theory and 0\t) is a trajectory 
on the unit 5*^. This operator is no longer unitary and so we cannot speak of an 
eigenvalue distribution on a circle. Instead, the expectation values (TrU'') determine 
an eigenvalue distribution on the complex plane which one can hope to use as an order 
parameter. For instance, it is this Polyakov-Maldacena loop which detects the strong 
coupling deconfinement transition. 

At large A^ and 't Hooft coupling A the AdS/CFT dictionary gives P^ITHITK] 

where S\k-wmding is the action of a classical string configuration in the dual bulk geom- 
etry that winds the thermal circle k times at the conformal boundary. Most previous 
work, both for temporal and spacelike Wilson loops, has focused on the single winding 
case k = 1. There are serious technical issues to overcome for multiply wound strings. 

A discussion of computing {TrU'') at strong coupling for spatial Wilson loops was 
given in |l5j. The approach is largely qualitative and is essentially an extrapolation 
to four dimensions of intuition that was developed in understanding two dimensional 
QCD as a simple string theory [ini HH d] ■ The classical string configuration is a disc 
worldsheet instanton embedded in the dual geometry as a fold cover of a minimal 
surface. The k winding requires the insertion of A; — 1 Z2 branch cuts running from 
interior points in the worldsheet to the boundary. This suggests the following expression 
to leading order in a' and string loop corrections jT^ 

^(Trf/'^) ~ (_i)'=-iA^-ie-'=^l'==i , (4) 

where the factor of (—1)'^^^ is due to the boundary conditions of fermions around the 
loop and A^~^ is a degeneracy factor due to the different possible locations of the k — 1 
branch points. 
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Despite the intuitive appeal of (jU it is not clear what A actually is. It could be 
the renormalised area of the infinite worldsheet or for instance, as suggested in (K], it 
could be the area of the worldsheet at the minimum of the background gravitational 
potential. Doing the calculation rigorously requires computing the measure on the 
moduli space of appropriate string instantons in the dual background. This does not 
appear to be technically feasible at the moment. Furthermore one might worry about 
the fact that the worldsheets have conical singularities at the branch cuts. 

In a beautiful paper, Drukker and Fiol have shown how to sidestep these issues by 
using a single probe D3 brane carrying electric worldvolume flux rather than fundamen- 
tal strings jTHI. In this approach there is no moduli space and worldvolume curvatures 
are small everywhere. They computed the value of k winding supersymmetric circular 
Wilson loops and gave a nontrivial match with a dual matrix model computation [20] . 
Even as oo, A ^ oo their result is significantly different from (0]) if k = k\fX/AN 
is kept fixed in this limit. In the following sections we will adapt their method to 
compute k winding Polyakov-Maldacena loops in the high temperature phase of the 
theory. We will find however that for the Polyakov loop case, nontrivial corrections 
come from D5 rather than D3 branes and are expressed in terms oi k' = k/N . 

Finally, we note that there is a perhaps underdeveloped parallel between this critical 
string theory story and two dimensional string theory on a circle. In that context the 
'Hawking- Page' transition is seen as a Kosterlitz-Thouless transition in which vortices 
condense on the string worldsheet pij. This condensation precisely corresponds to the 
insertion of branch cuts running to the boundary that we described above. Some of 
the vortex condensates have been computed explicitly using the dual large N matrix 
quantum mechanics on a circle P^. 

In section 2 we consider potentially dual D3 branes in the black hole geometry. 
We discuss possible boundary conditions and evaluation of the action on the solutions. 
We find that the only solution to the D3 brane equations of motion is a solution in 
which the D3 brane is collapsed on the cigar submanifold of the black hole. These 
configurations do not lead to the expected nontrivial k dependence. In section 3 we 
move on to consider D5 branes. Here we will find noncoUapsed solutions that capture 
the k dependence of the dual Polyakov loop. The two main results of this section are 
firstly the compution of the actions with /t' held fixed as ^ oo and secondly the 
observation that there are interesting A^~^/^ corrections to the large N limit with k 
fixed. Section 4 is a summary and a discussion of possible future directions. 
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2 (No) D3 brane solutions 



In this section we search for probe D3 brane duals to the multiply wound Polyakov 
loop. We will show that there are none, leading us on to a study of D5 branes in the 
following section. 



2.1 Equations of motion 

The dual geometry for Af = 4 SYM theory on x in the deconfined phase is the 
Euclidean Schwarzschild black hole in AdS^^ x S^: 



ds^ = R 



dr 



f(r)dt^ + -— + (da^ + sin^ adnl) + d^l 



(5) 



where 



fir) = l- '^^\^/^' +r\ (6) 



Following we are looking for a probe D3 brane configuration in the dual back- 
ground with the appropriate symmetries and worldvolume flux to describe a multiply 
wound Polyakov loop in field theory. The action for the probe D3 brane is a sum of 
the Dirac-Born-Infeld and Wess-Zumino terms 

S = TD3j drd^ae-'' ^det {*g + 2na'F) - ig^Tos j 'C4 , (7) 

where the tension = As usual F is the worldvolume gauge field and 

is the pull back of the bulk Ramond-Ramond four form potential. We can take the 
relevant part of the potential to be 

C4 = sin^ adt Ada A volS^ , (8) 

9s 

whilst the dilaton is contant: e~* = 1. 

The probe brane configuration must have the same symmetries as the dual Polyakov 
loop: SO (3) X SO (2) C SO (A) x ^0(2), the isotropy group of a point times S^ in 
S^ X S^. If {a,9,(f)} are coordinates on the S^, then the required configuration lies 
in the directions t,a,6,(j),r, with a nontrivial dependence on only one worldvolume 
coordinate 

a = a{a) , r = r{a) . (9) 

The geometry of this embedding is the cigar of the Euclidean black hole times an 
5*^ C S*^. Thus a(cr) G [0, tt] determines the size of the noncollapsed 5*^ inside the 
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S^. We can think of the D3 brane as k fundamental strings blown up via the dielec- 
tric Emparan- Myers effect ^Hl I24j . The string charge is induced from the worldvol- 
ume gauge field which must have electric field strength with nonvanishing component 
Fra{(j). The coefficient of the induced B field coupling to the string worldsheet will be 
the momentum 5S/5Ft-„. This must be set equal to ik to be reinterpreted as the charge 
of k fundamental strings. The presence of i will translate into the field strength being 
imaginary as was found in |19J. Therefore we introduce the notation Ft-o-((t) = iF{a). 
Evaluated on the ansatz we have just described, the action becomes 



S 



2N 



drdar'^ sin^ 



a 




fdaV 47r2F(a) 



A 



,da 
da 



(10) 



where we used = \a'^. 

The equation of motion for the gauge field gives the total string charge k E'L 

AN 27rr2 sin^ aF 



SS 

^ = -Jf 



A 



The equation of motion for a can be written 



:iii 



2 • 2 

r sm a 



■ 2 2nF ^ dr 

r sm2a— — -7- + 4r— 
kX^I^ da 



.kX^/"^ da 
2nF d^ 



where we introduced 



The equation of motion for r is 
27rF 



2r sin a 



sin^ a- 



/tAi/2 

One can use ()11|) written as 



^da 
da 



_ kVX 

1 kX'/^ djr^f) fda\ 

2 2tiF dr \db J 



d 

da 



kA1/2 
2'kF 'da 



/ da\ ' 
\d^) 



fdr\ 



(27rF)- 

K^X 



4 • 4 \ 

r sm a) 



(12) 
(13) 

(14) 
(15) 



to eliminate F from (|12|) or (jl4j) . The equations depend on F only through the com- 
bination 

G = ^. (16) 



kA1/2 ' 

SO we will work in terms of this quantity in what follows. 

Reparametrisation invariance of the action suggests that we should be able to con- 
sistently set 

r = a or a = a , (17) 

and then solve for a{r) or r(a), respectively. Using the three equations of motion above 
it is straightforward to show that indeed we may do this. 
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2.2 Boundary terms and boundary conditions 



Once we have solved the equations of motion, we need to evaluate the action on the 
solution to obtain the expectation value of the k winding Polyakov-Maldacena loop 



As is discussed very clearly in [19 , in evaluating the action on the solution it is im- 
portant to add the boundary terms. These terms implement Legendre transformations 
that permit the solution to have the correct boundary conditions: fixed winding number 
k and fixed momentum in the r direction. 



For our configurations we need to add 



dr r- 



6S 



Using the action ^ this term becomes 



dr At 



6S 



(19) 



r=r+ 



^1 



bdy. 



2N 

IT 



dr 



r dr 
G da 



2N 



TC 



drdan G . 



(20) 



In the above expressions, and in this section in general, we have assumed for simplicity 
that the configuration reaches the horizon r_|_. However it is possible for the solution 
to close off at a finite value r^aui > by having a(rmin) = 0. We will consider this 
possibility in later sections. 

There is a further ambiguity in the boundary term following from the fact that 
large gauge transformations of the background four form potential can change the 
Wess-Zumino part of the action. However, it was found in |Tn] that using a minimal 
expression for the potential compatible with the symmetries of the problem, as we have 
done, together with no extra boundary term gives the correct answer for cases that can 
be matched. 

In order for a D3 brane configuration to contribute semiclassically to a dual Polyakov 
loop it must have finite action, including the bulk and boundary terms. Let us look at 
the large r behaviour of solutions. We have found two asymptotic behaviours for the 
fields a{r) and G{r) at large r that lead to a finite total action 



KG(r) 



1 AB 

1 + 77-7-^^ + 



2r^ 



3 r-^ 



TC 



1 + 65^ 
6 



+ 



(21) 
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where B is an undetermined constant, and 



9A^K^ 3A^{4:Ak + 3) 
A A{2Ak + 3) 



^G{r) = 1 + ^- h^, '--^ + 



a{r) = TT 3 + , 3 ' - + ■■■, (22) 

where A is a free constant. A dependence on the horizon size r+ enters at higher 
orders in these expansions. We have written the expansions where a — > vr as r — > cxo. 
There are very similar expansions for a ^ 0. The a ^ vr case is more relevant as 
configurations with increasing a have a lower action because of the Wess-Zumino term 
in (fTUI) . In fact this term causes the action to behave very differently in the cases 
a ^ TT from below and a — > from above. Whilst both ()21|) and (j^^ lead to finite 
action as a — vr, only the faster falloff (j22|) is admissible as a — » 0. 

The falloffs (PT|) and (j^^ are special one parameter subfamilies of the general 
behaviour of solutions near infinity, which depends on two constants of integration. 
The generic behaviour at large r is 

^ , , mir) — rm'ir) 

a r = TT ^ + ■ ■ ■ , 23 

r 

where mir) ^ oo as r — oo but mir)jr 0. The explicit leading order behaviour 
of mir) appears to be somewhat complicated. Numerically we find that mir) ~ logr 
often gives a good description, but is not exact. However, without solving the equations 
for mir) it is possible to show that the action of these solutions behaves for large r as 



^\soin. = — I — Z:73 ZZTTIV + ^ • (24) 



-2Nf3 dr m(r)2 
n J 2 mi^r) — rm'i^r) 

This divergence may be very precisely checked numerically. Therefore configurations 
with the generic behaviour have infinite action and do not contribute to the 
Polyakov loop computation. It is not that the action is becoming unbounded below 
but rather that these boundary conditions do not give rise to any normalisable states. 

The bulk action is in fact finite with the asymptotic behaviours fIT\ and with 
another finite contribution to the action coming from the combined boundary terms. 
The action evaluated on a solution with all boundary terms included may be written 



2Np f 



S\soin. = —^y-r+^jK^ + rXsiYi a 



+ 



r'^fir) ( da\^ a o da r dG 
— — r sin a— — h — 



G \dr J dr G"^ dr 



dr (25) 
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is the period of the time circle. The simplest solution is the collapsed solution which 
has a = or a = TT. In this case the action is given immediately as 



\ collapsed ^ 2yj- ' \ ) 

This is k times the action for a fundamental string instanton wrapping the cigar. This 
connection between collapsed D branes carrying electric fliix and fundamental strings 
has been consistently verified in a range of examples since )23j . 

At this point we should note that for a general solution the dependence of the 
Polyakov loop on N and A is significantly restricted to be of the form 

S\soln=NPs{K,T^), (28) 

for some function r+). As was emphasised in ^H], although the expression (j^H|l 
is leading order in 1/iV and A, it is exact in k = kVX/AN. This represents a partial 
resumation of bulk higher genus string corrections. Such k dependence could not 
have been seen in the heuristic arguments of J3] leading to (0)) that we reviewed in 
the introduction. The result of for the k winding circular spatial Wilson loop, 
Sk = —2N [niy/l + + sinh~^ k\ , reduces to the naive expectation of Sk = kSi only 
in the limit k, 0. 



The important lesson from [in] is therefore that the possibility of a D3 brane blowing 
up due to an Emparan-Myers effect corresponds to and captures the nontrivial A and 

corrections which one expects from the singular nature of the fundamental string 
picture. These corrections are contained in the potentially nontrivial dependence of 
llzali on K = kVX/AN. We will now see that in contrast to the spatial Wilson loop case, 
this potential is not realised for the Polyakov loop. Instead, we will need to consider D5 
branes in the following section and find that these do give corrections to the collapsed 
result in terms of k' = k/N. At large A these are subleading compared with the k 
corrections which appear to vanish. 



2.3 Solving the equations 

After extensive numerical investigation of the equations of motion ()12j). ()14|) and ()15p. 
it has become clear that the only solutions which run from infinity to the horizon are 
the collapsed solutions a = and a = tc, with action ()27|1 . In this section we will 
combine analytic and numerical arguments to understand this fact. 
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We begin with some analytic results that are possible at large values of k. We will 
see in this regime that noncollapsed solutions do not reach the horizon but rather turn 
around at some minimum radius r^, and then go back out to infinity. We will then review 
continuity theorems for the solutions as a function of k and the boundary conditions to 
show that this behaviour persists for a range of solutions away from the large k limit. 
Finally, we present numerical data indicating that none of the noncollapsed solutions 
reach the horizon or close off in the interior. 

We will jump between considering r{a) and a{r) depending on what is most con- 
venient. G{a) or G{r) can be completely eliminated from the problem so we will not 
discuss its behaviour explicitly. 



2.3.1 The solutions at large k 



The equations of motion simplify when k is taken to be large. One may find the general 
solutions that have asymptotic boundary conditions (j2H) and (j22j) . With the former 
boundary condition the solution is 



r(a) = K 



l + B{a - vr) 



sm a 



(29) 



which can of course be viewed as a transcendental equation for a{r). In the latter case 
we have 

3A dy 
a{r) = 71-— . . (30) 



The large k expansion is valid for both of these solutions provided the constants B and 
A are of order one and for the coordinate range r+ <C r -C k^. Both the expressions 
(I29|) and ()30|) agree excellently with numerics. 

If we think of these solutions as r(a), then we see that as the solution comes in 
towards the horizon from r — oo it turns around at some minimal radius r* ~ k ^ r+. 
The behaviour near the turnaround point is 



r(a) = + (a - a^,)^ + ■ 



(31) 



The values of the constants r=K,a* and C may be determined as functions of i? or A 
from the solutions (j^^ and (jHUj) . respectively. If A and B are order one, then the 
turnaround occurs within the regime of validity of the solution. We note however that 
(I30|) will only be valid for the branch of the solution going from a = vr to a = a^,. 

None of these noncollapsed solutions reach the horizon and close off. Instead they 
come in from infinity at a = vr and turn around and go back out to infinity at a = 0. 
The solutions are illustrated schematically in Figure 1. The interpretation of these 
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solutions is as potential contributors to the two point function of the Polyakov loop at 
antipodal points of the field theory 5''^, that is (trf/''''^(0)trf/'^(7r)). However all of these 
solutions have an infinite action due to their behaviour as a — > 0. It is likely that the 
dual configuration for this antipodal two point function will be the disconnected sum 
of the one point functions, corresponding to a total screening of external quarks. This 
phenomenon is known to occur for the singly wound case fIE\ I26| ITTj. 




Figure 1: The solution reaches down to a minimal radius r^, > and then goes back 
out to infinity. 

Thus, at large k at least and provided we don't take A or i? to be very small, 
the only solution which could contribute to the expectation value of the k winding 
Polyakov loop is the collapsed solution. We will now argue that this continues to be 
true for all values of k and all values of A and B. 



2.3.2 Continuity as a function of k and boundary conditions 

To study the existence and properties of solutions to ordinary differential equations, 
the key question to establish is where the equation satisfies the Lipshitz property. To 
do this we write the equation for a(r) as two first order equations 

da 

dr ^ ' 

^ = F{a,^,r), (32) 

where the explicit expression for F(a,7,r) is easily computable and not particularly 
illuminating. The Lipshitz property in this case boils down to requiring that for all 
t, «!, a2, 7i, 72 in some rectangle we have 

71, r) - 71, r)]2 < L ([«i - aa]' + [71 - 72]') (33) 

for some constant L. Where the Lipshitz property holds, solutions to the equations 
exist, are unique given initial values of a and 7, and depend continuously on initial 
conditions and on parameters in the equation f2E\ . 
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For our equation ()32p the Lipshitz property may be seen to hold everywhere except 
when r— >r_|_,r^ooor7^oo. Therefore, away from infinity and the horizon, the 
only way the solution can terminate at some finite is if the derivative a'{r) diverges 
at that point. We have seen that indeed this occurs in the large k solutions, as (p?T|l 
implies a ~ {r — r^y/"^. One may show that such square root behaviour is the only 
allowed behaviour as a'{r) — > oo. 

Continuity then implies that as we vary k and the constants A, B in the finite 
action boundary conditions ()21|) and ()22j). all the solutions coming in from infinity will 
turn around at some finite r,., so long as does not approach the horizon r_|_. As we 
approach r+ the continuity argument no longer holds. Another possible termination is 
that a;(rmin) = for some rmin > r+ which would allow us to truncate the solution at 
that point as the S"^ closes off there. 

Although continuity therefore gives us a window around the large k solutions in 
which the behaviour is qualitatively similar, it is not enough to cover all possible /t and 
initial data. To this end we now present some numerical results. 

2.3.3 The minimum radius as a function of boundary conditions 

Figure 2 shows how the angle a,, at which the turnaround occurs varies as the minimum 
radius r* comes close to the horizon r+. The plot is obtained by varying the boundary 
conditions at infinity. A solution will reach the horizon if r^, ^ r+. However, what we 
see happenning is that as — > then a* — ^ vr. That is, as the solution gets closer 
and closer to the horizon, it tends towards the collapsed solution. The behaviour is 
illustrated in Figure 2 for three values of the horizon size, but in fact occurs for all 
horizon sizes from thermal AdS, which has r+ = 0, to the large horizon limit r_|_ oo, 
which may be described using a planar horizon. 

In Figure 2 we have used solutions that have the falloff ()21|) as r ^ oo. We also 
performed the computation with the faster falloff (f^^ and found exactly the same 
behaviour. As a crosscheck of these results, we also integrated numerically out from 
the horizon. We found that the solutions always reach r — >■ oo with a non normalisable 
falloff. The behaviour at infinity tends towards the normalisable falloff as Q;(r+) — vr, 
in perfect agreement with the results we have just described for integrating in. 

Therefore, we seem to find strong numerical evidence that all noncollapsed solutions 
satisying normalisable boundary conditions as r — oo do not reach the horizon. We 
have also seen that they do not close off by reaching a = at some rmin. It follows that 
the only D3 brane configurations that could contribute to the dual multiply wound 
Polyakov loop expectation value are the collapsed solutions. 
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Figure 2: The value of a* at which the turnaround radius r^, is reached for three 
values of the horizon size r+. From bottom to top r+ = 0.2, 8, 30. All plots have k = 5 
and asymptotic behaviour ()2H) . We find that a^, ti a.s r^, —>■ r^. 

We also considered whether any of the noncollapsed solutions gave a finite action 
contribution to the two point function (trt/^^(0)trt/^(7r)). It seems that none of the 
solutions have normalisable falloffs at both a — vr and a — > 0. 

2.3.4 Implications of the D3 brane results 

In this section we have studied D3 brane embeddings into the Schwarzschild AdS 
background, searching for dual configurations to multiply wound Polyakov-Maldacena 
loops in the strongly coupled large thermal field theory. The impressive matching 
of [1^, for the case of multiply wound spatial Wilson loops, indicated that the blowing 
up of an S"^ in the D3 brane worldvolume corresponds to a departure from the naive 
expectation Sk = kSi for the action of a multiply wound string worldsheet. This 
imphed that for spatial Wilson loops -^(Trf/^) 7^ -^(Trt/)^, generically. We have 
found that this blow up does not occur in the case of the Polyakov-Maldacena loops at 
strong coupling. If this were the full story it would imply instead ^^{TyU^) = ■^{TiU)'' 
following directly from ()27|). 

As we have discussed above, the Polyakov-Maldacena loop differs from the Polyakov 
loop by the presence of scalar fields and is not unitary. A priori one might expect A^^ 
degrees of freedom and hence obtain an eigenvalue distribution on a plane rather than 
a circle. However, the j^{TTU^) can be used to define eigenvalue distributions on a 
circle. Writing z = re*^ and fixing the modulus r we can obtain a distribution on the 
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6 circle 

oo ^ 

Pr{e) = l + 2'^r^coske—{TiU^). (34) 
fe=i 

When j^{TtU^) = -^(Trt/)*^ there are three possibihties, depending on the value of 
a = r^{TTU). If a = 1 then the distribution is a delta function. If a < 1 then the 
distribution becomes 

1 - 

p(e) = . (35) 

' l + a2-2acos^ ^ ' 

This is a smooth ungapped distribution on the 9 circle. If a > 1 then the sum does not 

converge. Analytic continuation in a leads to the result but now the distribution 

has a singularity that is not integrable. 

It remains to be seen what the precise relation of these distributions is to the eigen- 
value distribution of f/, and whether they can be used as order parameters. Further- 
more, the result j^(IiU'^) = -^(Trf/)'^ is clearly not reliable. The initial motivation 
for studying D3 branes was precisely to avoid issues that arise in the multiply wound 
string picture. However, the fact that the D3 brane configurations are collapsed seems 
to bring us back to that picture. Although the collapsed D3 brane retains a finite 
tension which is exactly k times the fundamental string tension, it remains true that 
induced curvatures on the collapsed S"^ are large and not controlled within the validity 
of the Dirac-Born-Infeld action. 

If the corrections to the multiply wound string picture are necessarily captured 
by a dependence on k = k^/X/AN then the fact that we do not find these corrections 
would be sufficient to imply that the eigenvalue distributions discussed above are indeed 
possible at leading order in and A. However, we will see in the following section that 
by considering D5 brane probes instead of D3 branes one finds corrections in terms of 
K,' = k/N. Note that in the large A limit ^ k and so these corrections are subleading 
with respect to potential D3 brane corrections. 



3 D5 brane solutions 

In this section we search for probe D5 brane configurations in the Euclidean AdS 
Schwarzschild background that have the correct charge and symmetries to contribute 
to multiply wound Polyakov loops. The logic will be close to that of the previous 
section so we shall be briefer in our presentation. 
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3.1 Equations of motion 



The action for the probe D5 brane is again a sum of Dirac-Born-Infeld and Wess- 
Zumino terms 



(36) 



where Tr)5 — Ny/X/Sn'^R^. As we will be considering solutions that are blown up in 
the direction, the relevant part of the four form potential is now 



C4 



3(7 - tt) 



— sm 7 cos 7 — cos 7 sm 7 



(37) 



where 7 G [0,7r] is a polar coordinate on the S^. 

Wc need a configuration that has the symmetries of the dual Polyakov loop: the 
isotropy group 5*0(3) x 5*0(2) of a point times in 5*^ x S-^, as well as an 5*0(5) C 
5*0(6) of the R symmetry group. This is obtained by having the D6 branc wrap an 
5*^ in the 5*^ and wrapping the time circle, while remaining at a point in the horizon 
5"^. Thus again the only nontrivial dependence is in 7(0-) and r(cr). The electric 
worldvolume gauge field will again be imaginary, so as before we let F^(t(c) = iF{a). 
Evaluated on this ansatz, the action becomes 



5" = — / drda 
37r2 



dr 



da 



sinSWh^ +f{r){-r 



d-fV 47r2F(c7)2 



da 



where we used i?^ = Xa'^ and introduced 



(7) = sin^ 7 cos 7 H — cos 7 sin 7 



A 



3(7 - tt) 



D(7) 



2nF{a) 
AV2 



(38) 
(39) 



We will set r — a and look for solutions for 7(r). The equation of motion for F is 



_SS _2N 
5F~^ 



2nF 



sin^7 



Ai/2 



(a - 



+ ^(7) 



(40) 



where as before A; e Z is the induced fundamental string charge. The equation of 
motion for 7 is 



4 sin 7 



sm 7 cos 7 



+ 1 



27rF _ d /AV^ 
AV2 ~ dr\ 27rF 



37r«;' 



fir) 



dj 
dr 



(41) 
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where we introduced 

It will also be convenient to introduce 

G=|^. (43) 

Note that the equations of motion are invariant under k —>■ N — k and 7 ^ tt — 7, 
reflecting the fact that the Polyakov loop is sensitive to the N-ality of the source. 

These equations are closely related to the D5 brane conflguration that is dual to the 
baryon vertex [SnUSOlE]- In fact, the solutions we will present clarify the uncertainty in 
these papers concerning the interpretation of solutions with k < N units of worldvolume 
flux. The simplest solutions we flnd will be at a constant angle 70 on the 5*^, similar 
in some regards to the conflning string solutions described in [3^. The connection 
between the baryon vertex and conflning strings was realised particularly explicitly in 



3.2 Boundary terms and action 



The boundary terms to be added to the action are similar to the D3 brane case. The 
action now depends on the derivative of one of the coordinates on the S^, 7. Therefore 
we need to include an extra boundary term to impose Neumann boundary conditions 
in this direction. A clear discussion of these issues can be found in [19j. The full 
boundary term we need to add can be written as 



^1 



bdy. 



-D(7) 



r + (7-7r)/(r)^ 



' 37r2 



drdr G . 

r+ 2 



(44) 



As with the D3 branes, there are collapsed solutions with 7 = vr. As in the D3 
brane case, these solutions take us back to the multiply wound string picture and are 
not reliable. These are again seen to have action 



^1 



collapsed 



\^kl3r+ 
2^ 



(45) 



which is k times the action for a fundamental string instanton wrapping the cigar. For 
a general solution we will have 



^1 



soln. 



(46) 
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for some function s(/€',r+). The D5 brane probes are sensitive to different corrections 
to the dual Polyakov loop than the D3 branes, depending on k' rather than k. Note 
that <^ /€ in the large A limit. 

As r cxD, the following falloff is allowed by the equations of motion and leads to 
a finite action configuration: 



7(r) 



TT 



A A(2 
r 



1 2A^ 1 1 
+ 7^ + 



(47) 



with A an arbitrary constant. This is the falloff considered by [SnilSni and, if the time 
circle were not thermal, would give an asymptotically supersymmetric configuration. 

There is another type of finite action solution which has qualitatively different 
behaviour. These are solutions in which the angle 7 tends to a constant 70, determined 
by k', as r — > 00. In particular, these include configurations in which 7 = 70 is constant 
everywhere. 

In contrast to the D3 brane case, numerical investigation of the equations reveals 
that integrating inwards with these boundary conditions leads to solutions that reach 
the horizon or close off at some radius where 7(rmin) = 0. Some possibilities are 
illustrated in figure 3. The case where 7(rmin) = is was considered by [211 ED]- In 
this case the solution is in the same homology class as the which therefore induces 
N units of fiux on the worldvolume. These solutions are only possible with k = N and 
correspond to baryon vertices in field theory. We are more interested in the solutions 
which reach the horizon as these can have arbitrary k and can therefore contribute to 
multiply wound Polyakov loops. These solutions can only exist because of the presence 
of a horizon, which corresponds to being in the deconfined phase in which nonsinglets 
can be screened. 









/ (3) 



Figure 3: Three types of D5 brane solutions. In (1) the solution runs from 7 = tt 
to 7 = 7+ at the horizon. In (2) the solution runs from 7 = tt to 7 = 0. In (3) the 
solution is at a constant 7 = 70. 

We postpone a systematic study of solutions to the D5 brane equations, and their 
field theory interpretation, for future work. Instead, in the following two subsections 
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we study firstly the simplest solutions, which are constant 7 = 70, and secondly we 
look at the k' ^ limit which is the most relevant for making direct contact with dual 
field theory computations. 



3.3 Constant solutions 

The configuration 7 = 70 is a solution to the equations of motion if 70 satisfies 

n^K — 1) = sin7o cos 70 — 70 • (48) 

This remarkable solution is possible because the metric function /(r) only appears in 
the equations of motion multiplied by dj/dr, and so drops out when 7 is constant. 
There is a unique value of 70 G [0, vr] for each value of k' G [0, 1]. 

The action, with boundary terms included, evaluated on these solutions is finite 
and may be written 

= ^2 — ^0 • (49) 

The action (j49|) is exact if k' is held fixed in the N —>■ 00 limit. This is not usually what 
one does in field theory, where the large limit is taken prior to computing Wilson 
loop observables, with the notable exception of the supersymmetric circular Wilson 
loop in A/" = 4 SYM theory which allows subleading corrections to be computed 
using a matrix model fK]\ We will consider the more usual N ^ 00 limit in 
the following subsection. It is also not clear that these constant solutions are the 
correct duals for -^(trt/'^), or any other higher loop in a representation of N-ality k, as 
opposed to the non constant solutions discussed around figure 3. However, given that 
the constant solutions may be characterised so explicitly, let us consider the eigenvalue 
distribution that would follow from using P9|l to compute -^(trf/^). As we note below, 
recent developments disfavour this interpretation. However, given the unexpectedly 
interesting result shown in figure 4, it seems possible that this computation will play a 
role in a more fully developed understanding. 

We will compute the eigenvalues numerically for a fixed, relatively large value of 
A^. If we were to assume that the result ()49|) gave us the traces 

^(Trf/'^) = e^^^'^'^'^oW, (50) 

then we may determine the eigenvalues by solving a system of equations. Here we 
introduced 

V = ^^. (51) 
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which is closely related to the action for a single fundamental string. We should take 
?7 ^ 1 corresponding to large 't Hooft coupling. At high temperatures jSr^ ~ 
Unfortunately, this makes the system of equations which determine the eigenvalues 
numerically rather intractible because the presence of Nt] in the exponent means that 
some equations have extremely large entries and others have entries of order one. 
We will not solve this problem here. Instead, let us take 77 ~ <^ 1. This of 
course takes us completely outside the regime of validity of the derivation of (jlUj) and 
leads to the important question of how strongly the eigenvalue distribution depends 
on 1]. Nonetheless, one obtains an interesting result. Figure 4 shows the eigenvalue 
distribution following from (j^Uj) with = 35. The eigenvalues lie on a circle in 
the complex plane, fully to the right of the imaginary axis. They are continuously 
distributed, with no gaps. 

2 [ . . . . . . . , 

1.5 

0.5 

ImA ' . 

-0.5 • .■ 

-1.5 

-0.5 0.5 I 1.5 2 2.5 3 
Re A 

Figure 4: Eigenvalue distribution of the Polyakov-Maldacena loop on the complex 
plane. With = 35 and N^t] = 1. 

3.4 Solutions for small k' 

This subsection looks at the limit — 0. This is the limit that makes the most 
immediate contact with standard field theory computations, where the strict large A^ 
limit is taken prior to calculating multiply wound loops. Hence k rather than k' is kept 
fixed as N ^ oo. 

To uncover the behaviour of the solution as — >• it is necessary to write 

7 = vr - (3^/2)^/^7. (52) 
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Then to leading order as k' ^ we have that 7 must satisfy the equation 



d 
dr 



d'j 
dr 



At (f - 1) 



(53) 



This is a polynomial nonlinear equation, in contrast to the general equations of motion 
we have been considering. The change of variables 7 = z^^^ and dr = f{r)dx takes the 
equation to a form that narrowly misses having the Painleve property. A very inter- 
esting result comes from computing the action on these configurations. The resulting 
action may be written 



S 



VXf3r+k 1 yA/5fc5/3 



27r 



+ 



iV2/3 307r2 



dr 



15f - 8f - 5/(r) ( ^ 



+ 



(54) 



The remaining higher order terms are not reliable in the strict large limit without 
including corrections to the background and the D5 brane action and so on. The 
first correction to the collapsed result appears to be reliable because A^"^/^ ^ N"'^. 
Numerical integration of (j53|) shows that there are solutions with finite action which 
run from infinity to the horizon. One further has the constant solution 7 = 1. If 
any of these provide the dual for multiply wound Polyakov loops or loops in higher 
representations, then within the validity of our computation we may write 



l + B 



03 



+ 



(55) 



where r] is defined in fl51|) and B is the integral given in In we have assumed 
a 1/r falloff which implies that boundary terms don't contribute. In the case of the 
constant solution 7 = 1 one has to add the boundary terms which result in a finite 
total action, as in the previous subsection (jl^ . 

The appearance of the non analyticity A^^^/^ suggests that the simple large A^ 
expansion is breaking down. In fact, (|55|) is incredibly tantalising because in the 
vicinity of a Gross- Witten transition ^U] , the large A^ expansion does break down 
and corrections need to be computed using a double scaling limit jHSl EEl EZj • Near 
the transition point, corrections are given precisely in terms of A^"^/^! Note that the 
nontrivial —2/3 scaling does not arise in the k ^ limit of the spatial Wilson loop 

m- 

We should comment on the validity of the solutions. There are two constraints: 
the requirement of negligible backreaction on the geometry and the requirement that 
worldvolume fields vary slowly |3H]- Adapted to the present situation, the former of 
these requires that k' <^ N and the latter requires k'^^^ ^ l/A^^^. The first is trivially 
satisfied while the second is certainly satisfied for the blown up solutions with finite k' 
and can be satisfied as ^ by taking A sufficiently large. 
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Note Added: Soon after this paper appeared on the arXiv two closely related 
preprints, ||J and were posted. In the zero temperature supersymmetric spa- 
tial Wilson loop context, these papers have clarified the interpretation of D5 brane 
solutions. It has been argued that these solutions compute traces of Wilson loops in 
the k-th antisymmetric representations of SU{N). Various interesting questions have 
been thrown up and we hope to address these in the near future. 

4 Summary and discussion 

In this paper we have studied the expectation value of Polyakov loops winding the 
thermal circle k times in A/" = 4 SYM theory at strong coupling. Inspired by the success 
of in computing multiply wound spatial Wilson loops, we began by searching for 
dual probe D3 brane configurations. We have shown that there are no appropriate 
solutions. Instead it appears that the correct dual description of higher Polyakov loops 
are probe D5 branes. The required configurations are, unsuprisingly, similar to those 
dual to the baryon vertex The blowing up of spheres in the brane worldvolume 
modifies the naive result jfiTiU'') = -^{TrU}^. In the D3 brane case the blowup is 
mediated by k = ky/X/AN while D5 brane blowups are determined hy k' = k/N. Only 
the latter corrections seem to be present for the Polyakov loop. 

There are various different probe D5 brane solutions of interest. We concentrated on 
two classes of solutions. Firstly, we studied solutions in which an S"^ in the D5 brane 
worldvolume is blown up at a constant size. The explicit nature of these solutions 
allowed us to present a numerical eigenvalue distribution. This distribution however 
had k' rather than k held fixed in the large limit and may not be directly comparable 
to the eigenvalue distributions computed in field theory. 

Secondly, we looked at solutions in the limit of small k,', or equivalently, N oo 
with k held fixed. We found that the leading correction to the naive multiply wound 
string result is of order N~'^/'^. This is a very curious result as N~'^^^ is the scaling one 
finds in the vicinity of a Gross- Witten phase transition. 

Thus while we are not yet in a position to make hard statements about the behaviour 
of the Polyakov loop eigenvalue distribution at strong coupling, and hence the phase 
structure of the finite temperature A/" = 4 theory, we have identified the required degrees 
of freedom on the bulk side of the duality and uncovered various hints of interesting 
behaviour. 

We are left with many further questions. One immediate open question is to sys- 
tematically understand the different possible D5 brane configurations and their field 
theory interpretations. In particular, one would like to know exactly which of the 
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configurations, if any, is dual to tlie multiply wound loop. 

The tantalising appearance of N~'^/'^ corrections suggested a connection with the 
double scaling limit near a Gross- Witten phase transition. A computation that could 
firm up this statement would be the calculation of (Trf/^), and other higher represen- 
tation loops, in the double scaling limit. One would be looking for corrections of the 
form A;5/3/7v2/3. 

It is worth bearing in mind that the Polyakov loop being a composite operator, its 
expectation value suffers ultraviolet divergences. The divergences in Wilson loops have 
been well understood ([22] and references therein) and result in an overall multiplicative 
constant whose origin is the infinite additive mass renormalisation of a static test quark 
in the theory. 

Finally, another possibly interesting direction for further investigation would be to 
compute the traces of the Maldacena- Polyakov loop in a partially resummed field theory 
model. Although this is not a controlled approximation for non supersymmetric loops, 
it has recently been seen to capture important physics away from the perturbative 
regime in the 't Hooft coupling |^ . 
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